Greenberger-Horne-Zeilinger generation protocol for superconducting transmon 

qubits capacitively coupled to a quantum bus 
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We propose a circuit quantum electrodynamics (QED) realization of a protocol to generate a 
Greenberger-Horne-Zeilinger (GHZ) state for A'^ superconducting transmon qubits homogeneously 
coupled to a superconducting transmission line resonator in the dispersive limit. We derive an 
effective Hamiltonian with pairwise qubit exchange interactions of the XY type, g{XX + YY), 
that can be globally controlled. Starting from a separable initial state, these interactions allow to 
generate a multi-qubit GHZ state within a time tcHZ ~ ■ We discuss how to probe the non-local 
nature and the genuine A''-partite entanglement of the generated state. Finally, we investigate the 
stability of the proposed scheme to inhomogeneities in the physical parameters. 
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I. INTRODUCTION 



Entangled quantum states are one of the essential re- 
sources for quantum information processing. They are 
necessary for the realization of quantum communication 
and the most important computational tasks. Many ef- 
forts have been devoted to the elaboration of physical 
systems enabling the generation and the control of such 
states. In particular, different types of superconducting 
qubits are promising candidates to solve this problem. 
Until recently limited to two qubitsi"— , efforts to entangle 
superconducting qubits have lately reached a new mile- 
stone with the experimental demonstration of three-qubit 
entanglement^i^. 

In the present paper we consider transmon qubit9^>2i in 
a circuit quantum electrodynamics architecture^— and 
present a way to generate GHZ statesi^, i.e., maximally 
entangled states. Although the mathematical descrip- 
tion of multipartite entanglement for more than three 
qubits is still debatedi^"— , GHZ states remain paradig- 
matic entangled states which are, in particular, useful 
for fault-tolerant quantum computing or quantum se- 
cret sharingi^. So far, many different protocols have 
been proposed to generate such states in circuit QED 
setups^"—. Some of them are of probabilistic nature, 
i.e., if a measurement on the iV-qubit system has a spe- 
cific result, the system is known to be in a GHZ state af- 
ter the measurement^^"— . In Ref.j^, a M0lmer-S0rensen 
type^ one-step scheme to generate GHZ states both for 
superconducting flux qubits and charge qubits was pro- 
posed. The procedure is independent of the initial state 
of the quantum bus and works in the presence of mul- 
tiple low-excitation modes. However, higher excitation 
modes of the quantum bus will introduce inhomogeneity 
because of the shorter wavelengths of the higher modes 
and decrease the GHZ fidelity. Moreover, uncontrolled 
dissipation might be coupled through the higher excita- 
tion modes and induce extra noise. It would be ideal 
to devise a GHZ generation scheme that, while keeping 
the one-step, deterministic nature, would involve only a 



single mode of the quantum bus mediating the qubit in- 
teraction. 



For this purpose, in the present paper, we consider 
N superconducting transmon qubits homogeneously cou- 
pled to a superconducting transmission line resonator in 
the dispersive limit, i.e., the architecture realized in a 
number of experimenta^'^iii'S^"— . We show that the sys- 
tem is characterized by effective qubit exchange interac- 
tions of XY type that can be globally controlled. Start- 
ing from a separable initial state, these interactions allow 
to generate a GHZ state in a deterministic one-step pro- 
cedure. We discuss how to probe the non-local nature 
and the genuine A'^-partite entanglement of the gener- 
ated state and investigate the stability of the proposed 
scheme to inhomogeneities in the physical parameters. In 
contrast to Ref. |23| . the qubit-resonator interaction does 
not commute with the free Hamiltonian, and the qubit 
frequencies are tuned close to one resonator mode. The 
time evolution of the system is described by an effective 
Hamiltonian which allows a direct implementation of the 
M0lmer-S0rensen idea. Our scheme is the first one-step 
deterministic generation protocol of GHZ states which 
could be possibly implemented in the currently available 
circuit QED design. 



The paper is organized as follows: in Section |TT] we 
derive an effective Hamiltonian for N transmon qubits 
capacitively coupled to a superconducting transmission 
line resonator in the dispersive regime. In Section [m] we 
describe the protocol for generating GHZ states in our 
system. In Section IIVI we discuss ways to confirm the 
A^-partite nature of the entanglement in the generated 
states, and in Section |V] we study the effects of non-ideal 
physical parameters like inhomogeneities in the qubit- 
resonator coupling constants. 
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II. FULLY CONNECTED NETWORK OF 
TRANSMON QUBITS IN THE DISPERSIVE 
LIMIT 

We propose a solid-state implementation, based on an 
architecture of superconducting transmon qubits capaci- 
tively coupled to a quantum bus and derive an effective 
Hamiltonian for the system, which exhibits the appropri- 
ate XY exchange interaction. 

Transmon qubits consist of a superconducting island 
connected to a superconducting electrode through a 
Josephson tunnel junction with capacitance Cj and an 
extra shunting capacitance Cb- A gate voltage Vg is ap- 
plied to the island via a gate capacitance Cg, allowing 
to tune the dimensionless gate charge Ug — CgVg/{2e). 
The system is characterized by the charging energy Ec — 

I (2Cs), where Cs = + Cj + Cb is the total capac- 
itance of the island, and the Josephson energy Ej of the 
tunnel junction. 

Such Josephson junction based qubits behave effec- 
tively as quantum two-level systems in different regimes, 
categorized by the ratio Ej/Ec- We will focus on the 
so-called transmon regime, when Ej/Eq ^ 50 — 100. In 
this case the Hamiltonian of a single transmon qubit 
can be expressed as 

n^h^AEcifi-ngf -Ejcosip. (1) 

In the following we assume that the Josephson junc- 
tions form a dc- SQUID i.e., Ej is tunable by an external 
applied magnetic flux <i>oxt allowing to control indepen- 
dently each qubit. In this case Cs = Cg + 2Cj^' + Cb 

and Ej = 2£;j cos{7r$ext/$o) with C^' and Ej the ca- 
pacitance and the Josephson energy of a single junction. 

If a qubit is capacitively coupled to a superconducting 
transmission line resonator, Cg is now the capacitance be- 
tween the superconducting island and the resonator. In 
that particular situation the gate voltage involves a de- 
part and an extra term depending on the state of the res- 
onator, Vg — Vg'^ + V{x). Therefore the interaction with 
the resonator appears via the gate charge Ug, which im- 
plicitly includes the voltage V{x). Transmon qubits are 
more robust to l//-noise than charge qubits due to their 
exponentially suppressed charge dispersion*. However, 
we assume that the gate of each qubit can be controlled 
separately by microwave pulses in order to perform single 
qubit quantum gates. 

For simplicity we consider the qubits to be coupled to 
a single mode of the resonator. This is a reasonable as- 
sumption if the qubits are nearly resonant with only one 
mode. Since higher modes have frequencies which are 
multiples of the fundamental frequency, we can tune the 
qubit transition frequencies such that the detuning with 
respect to one mode of the resonator is one order of mag- 
nitude smaller than the detuning to all the other modes. 
Under these conditions we can realize the dispersive limit 
for a single mode of the resonator and neglect the influ- 
ence of higher modes, as is the case in experiments using 
one transmon qubitS^. 
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FIG. 1: (Color online) Sketch of a possible coplanar geometry 
for the proposed device with TV = 4 qubits. Qubits (small blue 
squares) are placed around the maxima of the electrical field 
amplitude (red line), i.e. near the center and the ends of the 
(quasi-)one-dimensional resonator (gray strip). The second 
mode of the electrical field (red arrows) mediates the qubit- 
qubit interaction. Input and output ports of the resonator 
are drawn in black. 



For instance the qubits could be mainly coupled to the 
second mode if they are placed near the ends or the center 
of the resonator, that is the positions where the electrical 
field amplitude is maximal. Such a possible geometry is 
sketched in Fig. [1] Following the procedure of canonical 
quantization of a (quasi-) one-dimensional superconduct- 
ing resonator—, the voltage across the resonator is given 

by 

n^) = y^cos(^|^)(a + at). (2) 

The length of the r esona tor is and its resonance fre- 
quency ujr — 27r/ -y/Lpc depends on its capacity c and 
inductance I per unit length. The position along the 
resonator is denoted by a; S [— Lo/2, Lo/2] and a(a^) rep- 
resent bosonic annihilation (creation) field operators. 

Following Ref. 0, the system, composed of the res- 
onator and N transmon qubits, can be described with a 
generalized Jaynes-Cummings Hamiltonian. This Hamil- 
tonian is expressed in the basis of transmon eigenstates 
|j)g, where the indices q label the transmon qubits, and 

for readability we define the operators H^*'' = \j)q{j\q, 
<y^'X = + llf = + 1)9019 and set = 1, 

N 

9=1 J 

(3) 

The qubits frequencies wj''* are presumed to be tun- 
able through external magnetic fields $cxti changing 
the effective Josephson energies of the qubits E'-f = 

2Sj''' cos(7r$|,^]./<i>o), and the coupling frequencies gj''' 
depend on the position of the qubits. Invoking the 
rotating-wave approximation, we have neglected rapidly 
oscillating terms. In the transmon regime, we can only 
keep transmon-resonator coupling coefficients for neigh- 
boring levels, since terms like \i)q{j\q for \i ~ j\ > 1 are 
comparatively small. Moreover in the large Ej/Eq limit 
asymptotic expression can be obtained for wj'^ and gj''* 
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in first order perturbation theorjii, 



Keeping terms up to second order in gj/ l^j, we obtain. 



^+2 



E 



(<?) 



^(6j2 + 6j + 3), 



5o V7Tlcos(^-^j, (4) 



.eC, 



(9) 



(9) 



1/4 



(9) 



(-3) 



Loc' 



This form of the couphng frequencies describes the 
situation shown in Fig. [TJ The amphtudes of these cou- 
phng coefficients ^j^-* can be assumed to be approxi- 
mately homogeneous if the positions Xq of the qubits 
satisfy |xq/Lo| — or 1/2, since the electrical field am- 
plitude decreases quadratically with the distance from its 
maxima and since the size of the qubits is typically much 
smaller than the resonator wavelength in realistic sys- 
tems. However even if close to the center or the ends of 
the resonator, the qubits should be placed sufficiently far 
apart to reduce direct inductive or capacitive qubit-qubit 
coupling. There are also other positions that the qubits 
can be placed in (e.g. nodes of higher modes). However, 
the homogeneity of the coupling constants is important 
in our approach and should be taken care of. 



In the so-called dispersive regime Igj'V^j*''! ^ 1' 
when transitions frequencies of the transmon qubits 

j'+i are detuned from the resonator frequency cj^, ex- 
citations of the resonator are virtual and the latter will 
rather act as a quantum bus mediating effective qubit- 
qubit interactions. The transition frequencies of the 



(9) 



(9) 



ujy,\ — Lu^f^ and 



,(9) 



^0 + 1 



^Yj+i — i^r- In this 



transmon qubits are defined as uj 
their respective detuning as A 

regime, eliminating the direct interaction between res- 
onator and transmon qubits to lowest order in g^^^ / 
we exhibit an effective qubit-qubit interaction. This 
can be seen by performing the canonical transformation 
e^'He~^ , where 



N 



9=1 



J>1 



a' a 



N 

E 

9=1 
N 

E 

9=1 



(9) 



J>1 



EE 

9#9' j,3' 



Am) ^ 

%' /„(9)„(9') , ^(9)„(9') 



(6) 

Here the ac-Stark shifts xf^ i the two-photon transition 



rates 77^"^^ and the effective qubit-qubit coupling coeffi- 
cient ^j^? are given by 

(q) ^ \9j I 

A (9) ' 



Aq) _ 1 3j 9j+l 



2 A (9) A (9) 
^J + 1 



,(9) 



(7) 



;;(99') 



^33' 



,(9) *(9') 
3 ^3' 



3 ' 3' 
2A^fA^P 



Two-photon transitions can be safely neglected since 
the parameters 77^'^'' are small in the dispersive regime. 
An effective Hamiltonian Hoff is now obtained by re- 
stricting our Hilbert space to the computational sub- 
space, that is the first two levels of each transmon qubit 
{|0),|1)}®^. In principle, the qubit-qubit interaction 
couples any states of the qubits with more than one exci- 
tations to states that do not belong to the computational 
subspace (e.g. for iV = 3, the state |110) or |111) will be 
coupled to 1 020) or |021)). However, the amplitudes for 
these mixing processes of computational states with such 
non-computational states are of order g"^ /{Ec^) and will 
be neglected^^. Under these conditions, 



,~,(9) 



9 

~(99') 



(8) 



9,9' 



N 



^ = EE 



^(9) 



H.c 



(5) 



where x^'^ = X^'^ - xl'Vs, ^r, = Hi - Hq, a, = a+ + 
Uy = — i(cr+ — (T_). The resonator and qubit frequen- 
cies get slightly rcnormalized, that is a; = a;^ — Xi'^V^ 
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and uJq'^i = ojq'-* — Xo''''- The Hamiltonian has the de- 
sired YX-form, provided that all qubits have identical 
parameters: that is all qubit and coupling frequencies 
are homogeneous, Wq^^ — Q, \go'''\ = 9, Aq''^ = A and 
= = g = sVA. Using Eq. Q we infer that 
^(q) = x = -gEc/iA-Ec) < g, where Ec = W01-W12 is 
the weak anharmonicity of the transmon qubits. As men- 
tioned earlier in Eq. Q the qubit transition frequencies 
can be made homogeneous by tuning the flux biases $oxt- 
From now on we assume the gj'^ are homogeneous. This 
is motivated by a promising new transmon architecture 
with tunable coupling that has been proposed recently^!. 
Inhomogeneous coupling constants will be discussed in 
Sec. El 

Previous GHZ state generation protocols based on ho- 



,19-21 



ne- 



modyne measurement of the transmission lim 
glected the effective exchange interaction because of the 
large differences in qubit frequencies. In our case, the 
qubit frequencies luq^^ are tuned to be identical using the 
flux biases, and this effective interaction plays a signifi- 
cant role in the generation of the GHZ state in a one-step 
procedure as shown below. 

If the qubit and coupling frequencies are homogeneous. 



the total spin operators 



(9) 



and their cor- 



responding Casimir operator j"^ ~ + jy + can be 
used to write the effective Hamiltonian in a very conve- 
nient form, 



(9) 



Evidently, Hcs is diagonal in the basis \J, Jz)<E)\n), where 
the states \J,Jz), describing the states of the N qubits, 
are the eigenstates of the operators and Jz with re- 
spective eigenvalues J{J + 1) and Jz. The states \n), de- 
scribing the state of the resonator, are eigenstates of a^a 
with eigenvalue n > 0. Since [H, J^] — 0, any eigenstates 
of will remain so under the action of this Hamiltonian. 
In the following, we will restrict ourselves to such states 
with J — N/2. For example states with all spins aligned 
in a particular direction belong to this type and are there- 
fore an appropriate choice for the initial state. Setting 
J — N/2 in what follows, we denote \ J = N/2, Jz) by | J^). 
The eigenstates of HcS are | Jz) (g) |n) with eigenvalues 
e(n, Jz) =ujn + g{N/2 + l)N/2 + {n + 2xn)Jz - gJ^ 



III. PROTOCOL FOR GENERATING GHZ 
STATES 

The effective Hamiltonian HcB allows to produce GHZ 
states by turning on the interaction for a definite dura- 
tion icHZ- It was shown in Refs. [22||23 that a Hamil- 
tonian of the type gJ^ will produce a GHZ state af- 
ter the time 7r/(2g), starting for instance from the state 
\ l)q. Implementation of such scheme in other qubit 
system has also been proposed2H.iS. The multi-qubit gate 



exp(z7r is sometimes referred to as M0lmer-S0rensen 
gate. 

We conveniently choose an initial state with all the 
qubits aligned in the same direction, that is the maximal 



superposition state j-^o) = ^„ 



+ \l)q)/V2. We as- 



sume that the qubits and the resonator are initially in a 
product state and the state of the resonator at t = is 
denoted pres, 



(10) 



Moreover I^/jq) — \ Jx — N/2) and can be expressed as a 
linear superposition of the states \Jz) (see Appendix,), 



1^0) = 



1 



N 

^ fc=0 



(^) \Jz = k-N/2). (11) 



We define p(t) as the density matrix evolving un- 
der the action of the time-evolution operator U{t) = 
exp(— iHcffi); where "Heff is the effective Hamiltonian 
Eq. ©, 



p{t)^U{t)p{t = 0)U\t). 



(12) 



We assumed that g/A 1 and therefore we have ne- 
glected the effect of the canonical transformation e'^ on 
the state vector. This turns out to be particularly useful, 
since U{t) is diagonal in the basis \n), we can describe di- 
rectly the time evolution of the reduced density matrix of 
the qubits pqh{t), obtained by tracing over the resonator 
state, 



Pqh{t) Tri.os[/3(i)] 

1 NT^ , , 



2JV 

n.k.k' 



"lPrcs|ri;V(fc)(fc,) 

\Jz = k-N/2){Jz = k'^N/2\, 



(*)-Vfc',„(t)) 



(13) 



where (pk,n{t) = k {fit + 2xtn + gt{N — fc)). 

The Greenberger-Horne-Zeilinger (GHZ) states we aim 
to produce are of the following form. 



|GHZ±) 




N 



\0)q- 



9=1 



V2 



(14) 

which are standard GHZ states up to single qubit rota- 
tions. These states can be expressed as a linear superpo- 
sition of the states \Jz) as well (see Appendix J: 



IGHZ^ 



N 

E 



.i-nk 



l±ie 

rW^/2 



{l)\Jz^k-N/2). (15) 



To see why a GHZ state is produced after some time 
^GHZ we consider the effects of either exp{igtJ^) or 
ex.p[igt{J^ — Jz)] (for N either even or odd) on the state 
I Jz = k — N/2). We establish that one of the two possi- 
ble GHZ states Eq. (|14p is produced when gt = 7r/2 by 
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noticing that 



l + ie'"(fe+f-i) 



V2 



N even , 



V2 



,-»f+»f [(fc-f )2-(fc-f )] 



, TV odd . 



Thus, a GHZ state is produced for every odd multiple 
of time icHZ • The shortest preparation time is icHZ = 
V(25) 

However the remaining term of the effective Hamil- 
tonian in Eq. the one which is proportional to J^, 
induces a collective rotation of the final state. The ro- 
tation angle depends again on N and the state of the 
resonator. The state Pqb(iGHz) is, 



Pqb(iGHz) 



,|pres|n) |GHZ(a„))(GHZ(a„)|. (16) 



Here, 



|GHZ(a)) = 



-ia.Jz 




and 2an/7r = (fl + 2nx)/g for N even. For N odd, 
2an/n = {n + 2nx)/g — 1, and the relative phase 
exp{iT:{N ~ l)/2) in Eq. ^ is changed to exp(i7riV/2). 

We notice that the produced states p(tGHz) is not ex- 
actly the state depicted in Eq. (IT4l) and therefore certain 
constraints on the angles a„ in Eq. (ITBl) are required to 
generate the proper state |GHZ^). At low temperature, 
only the ground state of the resonator is significantly pop- 
ulated and (0 1 Pros 1 0) ^ (njprosl'T-) for n > 1. Thus we can 
restrict our considerations to a„=o and this translates to 
some condition on the ratio fl/g. 

To illustrate this we consider the resonator to be ini- 
tially in its ground state pics = \n = 0){n = 0\. The state 
|GHZ^) is indeed produced at ^ghz, provided we can 
tune the frequencies D, and g such that 



- = 4m + 2 

g 



N, 



m G Z . 



(18) 



If the above condition cannot be satisfied, some correct- 
ing pulse expliS^Jz) can be applied to the final state 
Pqb(^GHz) to obtain a proper IGHZ"*") state. The appro- 
priate pulse length Sn depends on N and the ratio fl/g, 



ON 



N -2] mod4 



(19) 



Furthermore 6n ^ implies Eq. ^ 

If not only the ground state of the resonator is 
populated, higher photon numbers n produce ro- 
tated GHZ states, according to Eq. (IT51) . We no- 
tice that (GHZ(a„)|GHZ(a„+fe)) = cos^(fc7rx/(2.9)). 



which means that if a IGHZ"*") state is produced 
for excitation number n, a slightly rotated state 
exp{—i'KxJz/g)\Gii'Z'^) is produced for n + 1 (since x < 
g). Assuming some correcting pulse expliSj^Jz) has been 
applied, the reduced density matrix of the qubits pqb is a 
mixture of rotated GHZ states with classical probabilities 
depending only on the initial state of the resonator. 



Pqb(iGHz)f 



|GHZ+)(GHZ^ 



(n|Prcs|n)e-^-"t-^^ IGHZ+) (GHZ+ |e' 



(20) 



n>0 



We will now show that it is possible to choose realistic 
physical parameters in agreement with our assumptions. 
Transmon qubits have typical frequencies f2/27r around 
10 GHz and coplanar waveguide resonators (the quan- 
tum bus) can be realized with frequencies a;/27r of the 
order of 10 GHz with high quality factors^. Transmon- 
resonator coupling frequencies g/2TT around 200 MHz is 
a reasonable assumption. Detuning the qubits from the 
resonator such that g/A ~ 1/10 would lead to an effec- 
tive qubit-qubit coupling of the order oi g = g/lO and to 
preparation time iGHZ of approximately 12.5 ns. 



IV. MEASURING THE GENERATED GHZ 
STATES 

The question of detecting and probing the states gen- 
erated in our scheme naturally arises. For iV > 4, there 
is no unique way to quantify entanglement. We will focus 
on a measurement of the Bell-Mermin operator'^^ defined 
as 



B = 



2i 



N 



N 



.(?) 



= 2^-1 (|GHZ+)(GHZ+| - |GHZ")(GHZ"|) , 

(21) 

whose expectation value for A^-qubit quantum states is 
bounded by \{B)\ < 2^^^, and the extremal values 
±2^^^ are reached for the states |GHZ^). The maximal 
value predicted by local hidden-variable theory is V2^ 
(VW^) for N even (odd), leading to an exponentially 
increasing violation for the states |GHZ^) with A^, the 
number of qubits. Therefore, a measurement of the Bell- 
Mermin operator leading to a result greater than V 2^ 
(\/2^~i) indicates the non-local nature of the generated 
quantum states. 

Other bounds can be derived for this operator: e.g., 
any separable state satisfies |Tr(p'^i3)| < 1. A signifi- 
cant bound can also be derived if the state is m-separable, 
i.e. describes a system which is partitioned in m sub- 
systems that only share classical correlations. In other 
words, a pure state is called m-separable, for 1 < m < A^, 
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if it can be written as a product of m states, 



i/Pi 



(22) 



where the {Pi} describe a partition of the N qubits. 
Thus, a separable state in the traditional sense is A^- 
separable. A mixed m-separable state p™ is defined as a 
convex sum of pure m-separable states, which might be- 
long to different partitions'^. Such an m-separable state 
satisfies Ti{p'^B) < 2^"™. Thus, any measurement of 
B with outcome above 2^~^ indicates that the state is 
not even biseparable {2-separable) and demonstrates the 
existence of genuine iV-partite entanglement. 

The Bell-Mermin operator expectation value can in 
principle be obtained experimentally. This operator can 
be expressed as a sum of parity operators, and inferring 
its expectation value would require 2^~^ parity measure- 
ments. 



N 



1=1 odd p 



IN-l 



N 



-rPil) 



q'=N-l+l 



(23) 

For each term, I is the number of factors ay and 
stands for the sum over the (^) permutations p that give 
distinct products. The states IGHZ"*^) defined in Eq. 
are those that give exactly ±1 for each of the 2^"-'^ terms. 

There are therefore 2^~^ parity measurements to real- 
ize which is possible only if one is able to generate GHZ 
states with high accuracy in a repeated way. Following 
Ref. [13, these parity operators could be measured by dis- 
persive readout. Since the frequency of the resonator is 



ac-Stark shifted lu 



it is possible to access 



the value of the operator . The value of the parity op- 
erator ai^^ can then be unambiguously deduced from 



N 



|a(«)) = (-l)f-^^ 



(24) 



Hence, we can measure all the needed parities by rotating 
the operators Cy*^ appearing in Eq. ([^ to tri*^ using 
single-qubit rotations. 

By means of Eq. (jl3p . we can give an expression for 
the time evolution of the expectation value of the Bell- 
Mermin operator, {B{t)) — Tr [i?pqb(i)]. For this pur- 
pose we can express the matrix elements of B in the basis 
of the states \Jz), which diagonalizes the effective Hamil- 
tonian. 



N 



B = 



E bk,k' \Jz^k'-N/2){J,^k-N/2\, (25) 



where 



k,k'=0 



Ok,k' 



1 

2i 



(26) 
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FIG. 2: Behavior of the function Gjv(iGHZ 4- t) for different 
A'', assuming for simplicity that 5n = 0. 



Hence, {B{t)) can be expressed as a sum of oscillating 
functions G^, indexed by the photon number n, 



(27) 



The functions are Fourier series over a finite 



range of frequencies w^/j/ defined as 
k')[{k + k' - N)g-n-2nx], 



^k,k' 



N 



k,k'=0 



where 



(fc 



(28) 



(29) 



Equation p7|) shows that {B{t)) is characterized by 
many oscillations on timescales ^ ^ghz, since the w^j,/ 
are of the same order as 3> X- However, the en- 
velope indeed reaches its maximum at icHZ, provided 
that only the ground state of the resonator is significantly 
populated. These fast oscillations are the manifestation 
of local rotations of the qubits, Eqs. (jlGHlTp . We have 
seen that this issue can be solved equivalently in two dif- 
ferent ways and that the state jCHZ''') is indeed gener- 
ated after tcHZ , either by applying some correcting pulse 
exTp{iS]\fJz), defined in Eq. (I19p . or by tuning the frequen- 
cies fl and g to satisfy the condition Eq. (|I8p . Assuming 
for simplicity that Sn — 0, we have then 



G'^(tGHz) 



„2N 



TTX 
I 

25 



,2^ 



TTX 
I 

2g 



(30) 



The fast oscillations of {B{t)) around tcHZ become 
sharper as the number of qubits N increases, as shown in 
Fig. [21 In the simpler case Sn — 0, the behavior of G% 
around icHZ is given by 



G? 



ivfeHZ 



r) 



T « 



(31) 
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FIG. 3: Temperature dependence of the maximum max(B) of 
{B{t)), for t ~ toHZ (squares), normalized by 2^~^ . The sohd 
hue shows the theoretical bound tanh(ix)/(2r)) for a resonator 
initially in the thermal state Eq. (|33[) . Inset: relative devi- 
ation e — tmax/icHZ — 1 of the time tmax at which max(i3) 
is realized compared to the predicted time fcHZ ~ n/{2g). 
Here we considered A*' = 4 qubits and the parameters are 
n/uj = 1.105, g/uj = 0.0105 and thus g/A ~ 0.1. Local 
hidden-variable theory only allows values of (B) below the 
dashed line. For A'' = 4 this value also corresponds to the 
upper bound for biseparable states. 



and that also means that we need a higher precision, for 
larger N , in controlling either the protocol time tcHZ or 
the correcting pulse. 

Finally, the maximal value (-B(iGHz)) can reach de- 
pends only on the initial state of the resonator pres, pro- 
vided the above considerations have been taken into ac- 
count. Equations ([^0]) and (PT|) show that 



2i-^Tr 



^{n\prcs\n) 

n=Q 



2N '^X\ ■ 2N '^X 

COS n — sm n 

2gy V 2.9 



(32) 

For instance, we assume pros to be a thermal state char- 
acterized by a temperature T, 



1 



n) (n| 



(33) 



In this simple case, the outcome of the Bell-Mermin 
operator measurement (i?(tGHz)) should be at least 
2^-itanh(w/(2T)). 

A numerical evaluation of {B{t)), using the Jaynes- 
Cummings Hamiltonian Eq. ([3]), shows good agreement 
with our theoretical analysis. We consider the ideal case 
of homogeneous qubit and coupling frequencies and we 
choose frequencies satisfying Eq. (|18l) such that Sn — 0- 
We look for the maximal value of {B{t)) around icHZ, 
that is for \t — tcHzl < and for the time imax at 



which this maximal value is realized. 
= 4 qubits are shown in Fig. [31 



The results for 



V. INHOMOGENEOUS COUPLING 
FREQUENCIES 

To estimate whether our scheme is robust against small 
random deviations in the physical parameters, we con- 
sider small inhomogeneities in the coupling strengths 
gj'^''. This effect will be investigated numerically and, 
for this purpose we compute the real-time evolution of 
the Bell-Mermin operator, using the Jaynes-Cummings 
Hamiltonian Eq. ([3]), truncated to the two lowest lev- 
els of the transmon qubits. This should capture the 
main features of this effect, since in our effective de- 
scription of the system Eq. the third levels of the 
transmon qubits only affect the ac-Stark shifts x^'^-' and 
renormalize the resonator frequency. Assuming the qubit 
transition frequencies are still homogeneous ujjf^ = 17, 
the inhomogeneity of the coupling frequencies ^q''' pro- 
duces inhomogeneous qubit-qubit couplings coefficients 
A. 



The coupling constants g^f"^ are assumed to be Gaus- 
sian distributed with mean g and standard deviation Sg. 
The notation {gq} denotes a particular set of coupling 

frequencies g^'^K The real-time evolution of the Bell- 
Mermin operator for one set of coupling frequencies {gq} 
is denoted {B^g^j{t)). 

For a given number of random realizations {gq} {ur 
around 200) with fixed Sg, we first calculate the mean 
value, 



n _ ' ^ 



Then, the maximal value (i3(imax)) defined by 



(B(W)) = niax(B(t)) 



(34) 



(35) 



is found. Finally the variances, above and below 
the maximal mean value (i3(imax)), of the particu- 
lar set {(i?{g^}(imax))} are calculated. The variances 
are calculated separately above and below, because the 
{B^g^j{t„^ax)) are not Gaussian distributed. We also cal- 
culate the median among the (S{g^}(imax)) and notice 
that the distribution is strongly asymmetric. 

Results for A^ = 4 and Sg/g between to 20 % are 
shown in Fig. 2] The time at which the maximum is at- 
tained is generally in good agreement with the predicted 
value ^GHZ = 7''/(2g), as long as g/A is small. The value 
of (-B(imax)) remains close to the ideal one for Sg/g of 
the order of a few percents and thus we notice that our 
scheme can tolerate some inhomogeneity in the coupling 
constants. 
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FIG. 4: Effect of inhomogeneous coupling frequencies (Jq'' 
witli mean g and standard deviation Sg. We sfiow tfie depen- 
dence of the maximal mean value (-B(fmax)) of {B{9<j}(i)) on 



Sg/g for t ~ fcnz (squares), 
dard deviation of {B[g^y{tn 



The error bars show the stan- 
:)) above and below the mean 



value. The median of (tmax)) (circles) is clearly above 

the mean value. Local hidden-variable theory only allows val- 
ues of (B) below the dashed line. For = 4 this value also 
corresponds to the upper bound for biseparable states. 



VI. CONCLUSION 

To conclude, we have shown that it is possible to gen- 
erate multipartite Greenberger-Horne-Zeilinger states on 
a set of transmon qubits in a circuit QED setup in a 
one-step deterministic protocol. In the dispersive limit 
g <^ A, such a system behaves as a fully connected 
qubit network with exchange interactions proportional 
to g = The preparation time of the protocol 

is inversely proportional to g. The non-local nature of 
the generated state can be investigated using a Bell- 
Mermin inequality. Moreover, we have derived and ap- 
plied bounds on the expectation value of the Bell-Mermin 
operator as a detection criterion for genuine A^-partite 
entanglement. Finally we have shown that our scheme is 
robust against small inhomogeneities in the coupling fre- 
quencies. The implementation proposed here looks like 
a promising way to generate GHZ states, and hopefully 
can be experimentally realized in a circuit QED setup. 
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Appendix: Schwinger representation of total spin 
operators 

We present briefly the Schwinger representation'^^ of 

the total spin operators Jx,y,z — \^q'^x'}y,z- This 
turns out to be particularly useful for calculations in 
the subspace of J^-eigenstates with maximal eigenvalue 
T (t ~^ ^) ■^here N is the number of spins. From now 
on we set J~N/2 and denote the states \ J — N/2, Jx,y.z) 

by \Jx,y,z)- 

States like \Jz) are sometimes referred to as Dicke 
states^-, they form a complete basis of symmetric A'^- 
qubit states, i.e., states invariant under any permuta- 
tion of qubits. We use for each qubit the standard ba- 
sis {|0),|1)} with the convention a'f^ 
ai'^|0), = -|0)„ 



1), = \l)q and 



N/2} 



|l)p(fe)|0)p(fc4 



|0)p 



(N) 



(A.l) 

with < A; < iV and where the sum is taken over the 
il) ~ k\(N-k)\ iioncquivalent possible permutations p 
that give different product states. 

The operators Ji are defined by means of two inde- 
pendent bosonic operators a and 6, with commutation 
relations [a, — [b, = 1 and [a, b] = [a, b^ = 0, 



Jx = ^{b^a + ah), 
Jz = ^(bh-a^a), 



(A.2) 



fulfilling the SU(2) algebra [J/, J™] Jn ■ Eigen- 

states of Jz can be expressed as 

where \na = 0,nf, — 0) is the vacuum state of the operators 
a and b. Since the choice of the operators a and b is 
not unique, we can equivalently introduce the operators 
c = (a — b)/V2 and d = {a + b)/V2, leading to Jx = 
^{d^d— c^c) and 



\J,Jx 



:|na = 0, 716 = 0) . (A.4) 



y/{J+Jx)KJ~Jxy- 

We straightforwardly obtain the decomposition of the 
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states \J, Jx) in terms of | J, J2) and in particular 



total spins, the Dicke states can also be expressed as 



N 



|J, = ±iV/2) 



|0>,±|l)g 



q=l 



V2 



\J, = k-N/2) 



-|n, = 0,nb = 0) (A.5) 



(at±6t)^ 
^ ^ I — 



|0>, 



N-k 



(A.' 



2N/2 



k=0 



Defining the ladder operators J± — Jx ^ iJy of the 
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